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In this paper, we present the recent enhancement of the Open National Combustion 
Code (OpenNCC) and apply the OpenNCC to model a realistic combustor configuration 
(Energy Efficient Engine (E°)). First, we perform a series of validation tests for the newly- 
implemented advection upstream splitting method (AUSM) and the extended version of 
the AUSM-family schemes (AUSM*-up). Compared with the analytical/experimental data 
of the validation tests, we achieved good agreement. In the steady-state E® cold flow 
results using the Reynolds-averaged Navier-Stokes(RANS), we find a noticeable difference 
in the flow fields calculated by the two different numerical schemes, the standard Jameson- 
Schmidt-Turkel (JST) scheme and the AUSM scheme. The main differences are that the 
AUSM scheme is less numerical dissipative and it predicts much stronger reverse flow in 
the recirculation zone. This study indicates that two schemes could show different flame- 
holding predictions and overall flame structures. 


I. Introduction 


Designing high-pressure turbines (HPTs) for peak temperatures at the combustor exit decreases cycle 
efficiency since excess cooling air is required, while designing HPTs for the mean exit-temperature at the 
combustor exit may reduce turbine life since excessive metal temperatures may occur at local hot spots along 
turbine blades. Improved modeling of combustor-turbine interactions can provide strategies for conditioning 
the combustor exit flow or optimized cooling strategies for HPT, which provides increased engine cycle 
efficiency and improves turbine blade life and overall gas turbine durability. Future propulsion systems will 
be of increasingly higher bypass ratio from larger fans combined with much smaller cores, which increases the 
importance of understanding core engine component interactions, such as combustor-turbine interactions. 
This work represents the first steps in a longer-term effort to investigate and quantify the effects spatial and 
temporal temperature non-uniformity of the combustor flow field on the high-pressure turbine performance 
using high-fidelity CFD simulations. 

Although there have been significant improvements recently in the field of computational fluid dynamics 
(CFD), the current maturity of the CFD capability of modeling a gas turbine and the relevant physics 
involved is not fully satisfactory. This is due to the fact that a gas turbine combustor is a complex system 
coupled with multi-physical phenomena, i.e. atomization, transport, vaporization and combustion, all of 
which are closely connected to turbulence. Design of gas-turbine combustors typically uses RANS during 
initial design cycles and Large Eddy Simulation (LES) for final design evaluation. RANS provides averaged 
solutions with relatively low computational cost (several hours for modeling a typical gas combustor). On 
the other hand, LES can provide an improved description of unsteady features, including turbulence and 
mixing, in contrast to traditional RANS methods.! LES also provides a tool to predict unsteady phenomena, 
such as flame stability, combustion dynamics, and combustor-turbine interaction (hot-streaks), all which are 
of our ultimate interest. 
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(a) (b) 


Figure 1. (a) Five-cup E® sector test hardware!? and (b) one-cup CAD geometry. 


At the National Aeronautics and Space Administration (NASA) Glenn Research Center, we have been 
developing a high-fidelity three-dimensional unsteady flow solver (OpenNCC), which is a publicly releasable 
version of the National Combustion Code (NCC).? A dual time-stepping approach is combined with a 4-stage 
explicit Runge-Kutta time-integration scheme to achieve the second-order accuracy in time. OpenNCC is 
capable of using unstructured meshes and has a massively parallel computing capability, with which we 
are able to simulate different types of combustors (e.g., Lean Direct Injection (LDI) combustion) related 
to a variety of NASA missions. For more than two decades, the team has continuously updated the NCC 
(and OpenNCC) by adding more physical modules, such as a turbulence module (nonlinear k — € model’), 
combustion modules (reduced chemical kinetic, low dimensional manifold), a spray module,*® and a joint 
probability density function for species and enthalpy. There are two numerical schemes implemented in 
OpenNCC: the standard Jameson-Schmidt-Turkel (JST) scheme and Roe’s upwind scheme. Recently, we 
implemented the advection upstream splitting method (AUSM)’ and the extended version of the AUSM- 
family schemes (AUSMt,® AUSM*-up’). The uniqueness of this type of scheme is that it does not have 
any computationally expensive matrix operation nor any differentiation of fluxes causing any numerical 
errors, yet it maintains accuracy. Another important feature is that there are not modeling parameters that 
need to be adjusted for applications (for instance, in the JST scheme, the second-order and fourth-order 
artificially dissipation coefficients need to be tuned, and the choice of these parameters are closely related 
to the numerical stability and then the final results). We have performed a series of validation tests for the 
newly-implemented AUSM scheme, such as skin friction at the laminar/turbulence flat plates and heat flux 
of the highly-loaded turbine guide vane at the Von Karman Institute (VKI), and methane-air reacting flow 
in a non-swirling coaxial jet combustor at the United Technologies Research Center (UTRC).!° 

This paper represents a first step toward the ultimate objective of tightly-coupled unsteady simulations 
of combustor-turbine interactions for realistic combustor and turbine geometries. To advance our current 
capabilities, we start with the non-proprietary combustor and high-pressure turbine designs from the Energy 
Efficient Engine (E*) program!! (using the General Electric combustor and HPT designs). While only the 
E? combustor geometry is considered in this paper, we plan to eventually simulate the E? HPT geometry 
and the coupled E? combustor-HPT geometry (shown in Figs. 1 and 2, respectively). The E® was intensively 
investigated in the mid-1970s and 1980s and set a historically important milestone toward more fuel efficient 
jet engines that meet the emission requirements. The relevant technologies that emerged from this program 
have since been widely used in current commercial engines. 

The E? is the double-annular and compact combustor. The fuel injected from 30 fuel nozzles and the 
swirling air are efficiently mixed and combusted in a short distance. At a low power condition, only the outer 
dome is fueled, and a rich combustion zone is formed. At a high power condition, both domes are fueled 
and a large amount of airflow is introduced into the inner dome annulus, creating a very lean combustion 
zone. Large dilution holes at the center body and inner/outer liners are designed in such a way that the 
cold airflow immediately mixes with the hot product, suppressing the further NOx formation. A detailed 
explanation of E® can be found elsewhere.!” 
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Although this type of combustor is designed for efficient fuel/air mixing and a relatively uniform exit 
temperature profile, temporary severe non-uniformity of temperature (i.e., hot-streaks) at the combustor 
exit is inevitable even for the latest and most advanced combustors. These hot streaks can cause local hot 
spots on the blade surfaces of high-pressure turbine stages, which is a potential source of blade life reduction. 
In order to improve understanding of these effects, numerous experimental and numerical studies have been 
carried out. However, due to the aforementioned coupled multi-physical phenomena, the numerical modeling 
is still not mature enough for a full understanding of the combustor-turbine interactions (and formation of 
hot-streaks). Thus, our goal is to enhance the current model capability of OpenNCC for this purpose. 

The layout of the paper is as follows. In Section 2, we briefly explain the OpenNCC and the current 
capability and subsequently show a series of validation works for the newly implemented numerical models. 
In Section 3, we show some preliminary work applying the current code to the E? and investigate the 
three-dimensional flow feature inside the combustor. Our conclusions are summarized in Section 4. 


II. Numerical Models 


A brief explanation of OpenNCC is provided in this section. OpenNCC is the releasable version of the 
NCC, which has been continuously updated for more than two decades at NASA Glenn Research Center 
(GRC). It is a state-of-the-art code equipped with a comprehensive combustion model, a turbulence model 
(cubic non-linear k—e€ model with the wall function), and a numerical simulation and spray model (Lagrangian 
liquid phase model). Many previous papers and presentations (e.g.,2-3-'417) have shown that OpenNCC (as 
well as NCC) has facilitated the combustion CFD in the development /design of the combustion technology 
at NASA-GRC. However, as with any other CFD code, there is a need to continue updating the code 
through verification and validation (V&V). Our experimental combustion group at NASA-GRC provides the 
validation data to improve the OpenNCC (e.g.,1%-19). For the numerical and mesh capability, OpenNCC 
is designed for unstructured grids (i.e., any mix of three-dimensional elements: hexahedral and tetrahedral 
mesh), and massively parallel computing (with almost perfectly linear scalability is achieved for non-spray 
cases up to 4000 central processing units). In addition, the adaptive mesh refinement option is recently 
implemented.?° For low-speed flow, pre-conditioning is used for better convergence. The combination of 
these unique features gives OpenNCC an advantage in the investigation of a variety of combustor /injector 
applications. 


III. Validation Test 


As a part of the continuous improvement efforts of OpenNCC, we recently implemented the AUSM* 
and the other extended version of the AUSM-family schemes (AUSM*,® AUSM*-up’) and the turbulence 
transition model.?! The AUSM scheme is particularly attractive for its very simple formulation (i.e., com- 
putationally cheap), no ad hoc model parameters, and great accuracy. Also, this simple formulation easily 
facilities the addition of other variables. The key idea is to split the inviscid flux, F, into the convective F, 


Figure 2. CAD Geometry of high-pressure turbine stages of E? (based on geometry data‘), 
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Figure 3. Skin friction (top) and velocity profiles (bottom) for the laminar flat plate. 
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Figure 4. Predicted wall friction on the turbulent flat plate and the experimental data. 


and pressure P terms, 


F=F.+P (1) 
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where F, and P are (p, pu, phy) and (0,p,0) for the one-dimensional system. p, u h; and p are density, 
velocity, total enthalpy and pressure. Then, the numerical approximations of the cell interface Mach number, 
the sound speed and the pressure are evaluated by the right and left cell variables. The extension of the 
high-order scheme can be done based on the MUSCL scheme.?? Please see other references for a detailed 
explanation of AUSM (e.g.,7°). A series of validation of the newly implemented AUSM scheme of the 
OpenNCC has been performed. The turbulence transition model developed by Liou and Shih?! on the 
platform of two equation models is also implemented. The model is based on the experimental observation 
that the transition takes place in the highly intermittent flow, so they introduce the intermittency correction 
function, y which reflects the peak disturbance energy. The implementation of the model into the OpenNCC 
is fairly straightforward; however, finding the peak disturbance energy in the wall normal direction within 
the massively parallel computing system should be done with care. 


A. Laminar/Turbulence Flat Plate 


The first case is to simulate the laminar flow over 
the flat plate. The results are compared with the 
Blasius solutions. The single block, two-dimensional 
structured mesh is generated by Cubit, a software 
toolkit developed at Sandia National Laboratories. 
The grid density is 128 x 64. The inflow and out- 
flow condition is: Majn, = 0.2, Tin = 300 [K], and 
Pout = 101,352 [pa]. The wall is treated as an adia- 
batic and non-slip wall. The resulting skin friction, 
u and v velocity profiles at the exit are compared 
with the Blasius solution (see Figs. 3 (a)-(c)). Ex- 
cellent agreement is achieved. 

The second validation case is to predict the skin 
friction of the turbulent flow over the flat plate. The 
experimental and empirical fit (C+ = 0.0576/Re°? 
are available. Two single blocks, two-dimensional 

Pseudocolor structured meshes (grid densities are 256 x 128 and 

-0.15 peg ae 512 x 256) are generated by Cubit, and we refine the 

mesh to adequately resolve the flow in the vicinity 

of the wall. These two meshes are used to check 

the grid dependency of the result. There are two 

nee operating conditions to test: Tu = 3.3 % and Tu = 

=Oe20 Taney a 6.5 %. The mixing length is chosen so that pu¢/p = 

Min: 4,593e-06 10. Fig. 4 shows the predicted wall shear stress 

using the 256 x 128 grid. There is good agreement 

-0.020.000.020.040.060.080.10 in the data with respect of the magnitude as well as 
X-Axis the transition location. 
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Figure 5. Contour of Mach number of VKI vane for B. VKI Heat Flux 


MUR 45. . . ; . : : 
The third validation test is to investigate the flow 


around the highly loaded turbine guide vane. The 

data was acquired at the Von Karman Institute 

(VKI), and the detailed test configuration and operating conditions are provided by the reference.?° They 
measured the velocity distribution by means of static pressure tappings and the wall heat flux by platinum 
thin films. Experimental uncertainties are also provided, pressure: 40.5 % and heat flux: +5 %. Here, 
we pick three free stream conditions, MUR45, MUR129 and MUR235 (see Table. 1) Note the inflow total 
temperature is not provided for MUR45. Thus, we estimate 410 [K] based on other experimental conditions. 
Figure 5 shows the contour plot of the Mach number for MUR45 (laminar calculation). We can see that 
the flow reaches supersonic speed at the suction side and then gradually decelerates toward the trailing edge. 
The current mesh allows the wake region to extend well in the computational domain without crossing the 
periodic wall (bottom and top walls), which minimizes the numerical errors. We observe the expansion on 
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Table 1. List of experimental free stream condition.?% 


Test number # | protat (bar)  Thotai[K] Misout Reout Free stream turb. [%] 
MUR45 1.475 - 0.875 10° - 
MUR129 1.849 409.20 0.840 1.1352 x 106 0.8 
MUR235 1.828 413.3 0.927 1.1521 x 106 6.0 


the pressure side in the vicinity of the trailing edge due to the large curvature. Figures 6 (a) and (b) show 
the pressure and the isentropic Mach number (Mis) with the experimental data (black circle) for MUR45. 
There are reasonably good agreements obtained, except the trailing edge region where we need to refine the 
mesh further more. Figures 7 (a) and (b) show the wall heat flux for MUR129 and MUR235. The region 
with negative X/C indicates the pressure side. For MUR129, the result using the laminar model agrees well 
with the experimental data and the one from Gourdain et al.?4 For MUR235 (see Fig. 7 (b)), both our model 
and Gourdain et al. (RANS) are similar to each other and predict much less heat flux at the laminar region 
than the data. In fact, Gourdain shows that using the high fidelity model (i.e., LES) performs better. Note 
that MUR235 is the most challenging case in the sense that there are transitions from laminar to turbulence 
as well as weak shock. 


C. UTRC Non-Swirling Coaxial Jet Combustor 


The fourth validation test is to simulate the coaxial jet combustor experiment conducted by Spadaccini et 
al.!° This configuration is relatively simple, but similar to a gas turbine combustor in the sense that in both 
cases, the flame is held by the recirculation zone. Detailed measurements (species, temperature and velocity 
profiles at different axial locations) are available. 

Figures 8 (a)-(c) show the computational domain and meshes. All the grids are structured and the mesh 
count is about 370,000.?° The inflow condition of the fuel and air are summarized in Table. 2. Figure 9 shows 
the axial velocity (left) and temperature (right) contours using two different numerical schemes, the JST 
scheme (top row)?° and the AUSM scheme (bottom row). In both cases the k-e turbulent model is used, and 
the chemistry-turbulence interaction model is turned off (i.e., the reaction rates are evaluated by the mean 
temperature). Compared with those utilizing the JST scheme, the AUSM scheme predicts a much stronger 
reverse flow, and as a consequence, the hot product stays upstream and the fuel stream expands more in 
the radial direction. We observe that the cold fuel core predicted by the JST scheme remains downstream 
in the temperature profile. Figs. 10 (a)-(e) and Figs. 11 (a)-(d) show comparisons of the axial velocity 
and temperature between the predictions by the AUSM scheme (red solid line) and the experimental data 
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Figure 6. (a) Pressure and (b) Mach number profiles for MUR 45. 
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Figure 7. Heat flux at pressure side (X/C<0) and suction side (X/C>0) for MUR45 (a) and MUR235 (b). 
The circle and plus are from?? and.?4 
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Figure 8. Geomtry of the coaxial jet combustor and meshes.?® 


Oo 


Table 2. Experimental setup of UTRC non-swirling coaxial Jet combustor. 


gas Temperature [K] Mass flow rate [kg/s] Velocity [m/s] 
Air 750 0.137 29 
Fuel (CHy) 300 0.0072 0.9 


(triangle)'° at different axial locations. From Fig. 11 (a), we can see that the hot product stays upstream 
in a way similar to that which is as predicted by the AUSM scheme. Overall, reasonably good agreement 
has been achieved. 
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Figure 9. Contour of axial velocity (left) and temperature (right) using two different numerical schemes: JST 
(top)?5 and AUSM (bottom) 
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Figure 10. Axial velocity profile at five different axial locations: z/d = 0.52, 0.146, 0.187, 1.58 and 1.79. The 
solid red line and triangle are the numerical prediction and the experimental data.1° 
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Figure 11. Temperature profile at four different axial locations: z/d = 0.34, 0.60, 1.73 and 1.99. The solid red 
line and triangle are the numerical prediction and the experimental data.!° 


IV. E°® Combustor 


In this section, we apply the OpenNCC to simulate the E? combustor and show some preliminary results 
using the RANS. 


A. Numerical Setting and Boundary Condition 


For this calculation, thanks to the massively parallel computing capability of OpenNCC, we routinely use 
960 processors of Pleiades at NASA Advanced Supercomputing facility. The tetrahedral mesh is generated 
by Cubit. The total mesh count is about 9.5 million. (Note that we are aware that using tetrahedral mesh 
for an internal flow should be done with care since it is numerically dissipative and not good at capturing 
the boundary layer). In future works, we would like to turn on the adaptive mesh refinement option in 
OpenNCC.?° 

For illustration purposes of the enhanced OpenNCC capability, all results shown here are obtained by 
using the RANS with the low-Mach preconditioning. The turbulence model is the nonlinear k—¢ model,? and 
two types of numerical schemes are used, the JST scheme and the AUSM scheme. In the JST scheme, we set 
the second and forth artificial dissipation coefficients to be 0.02 and 0.06, respectively. The E* experimental 
setup is shown in Table. 3 . The sea level takeoff condition (SLTO) is the most severe condition during the 
engine operation cycle. The pressured air (P3 = 2.52 [atm]) is introduced from the left boundary (A1) in 
Fig. 12. The air goes around the diffuser splitter vane, and some portion of air (S1) is diffused into the strut 
cavity and routed through the hollow strut passage into the cooling circuit. In the reacting flow case, both 
the main and pilot domes should be equally fueled (F1 = F2 = 0.00182 [kg/s]). There are several cooling 
airflow (A2 - Al4) at the outer/inner liners and the center body surface. We treat them as source/suction 
terms (i.e, no grids to resolve the cooling holes), and the mass flow rates are summarized in Table. 4. 
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Figure 12. Cross section of the sector E* combustor and injectors. 


Table 3. E® experimental setup.!? 


P3 [atm] T3 [K] W3 [kg/s] Wtota [kg/s| f/a  Wfpitot/Whtotat Tuer [K] 
SLTO 2.52 720 0.26 0.00364 0.014 0.5 520 


B. Cold Flow Comparison 


In this subsection, we show the cold flow (the domes are not fueled) results using the JST scheme and the 
AUSM scheme. ‘The contours of gage pressure, u-velocity, v-velocity and w-velocity at the mid-plane are 
shown in Figs. 13 (a)-(h). Both models are able to qualitatively capture the important features of the 
internal airflow inside the combustor in a similar way, such as the central recirculation zone (CRZ), high 
pressure region behind the domes, swirling flows, strong cross airflow from the dilution holes, and cooling 
air from liners covering the combustor liner. The predicted pressure drop across the combustor liner is 5.1% 
(main) and 4.7% (pilot), which is consistent with the designed pressure drop, 5%. The CRZ, which is formed 
due to the breakdown of a highly-swirled flow, plays a critical role in holding the flame in this type of a 
swirl stabilized combustor and is often associated with the low pressure region in the vicinity of the fuel 
nozzle, the so-called precessing vortex core (PVC). The inherent unsteady feature of the PVC significantly 
impacts the combustion dynamics and emission.?°?” Even though the presence of the PVC can improve the 
mixing because of the creation of larger turbulent scales, it is important to suppress the PVC in order to 
reduce the combustion dynamics, which can lead to catastrophic failure of the combustor. In the contours 
of u-velocity and v-velocity (Figs. 13 (c) - (f)), it is found that there is a large CRZ (blue region) located 
between the dome and the dilution flows, where the flame is supposed to be held in a reacting flow case. 
The wake region behind the dilution airflow coming from the center body can be a potential region where 
some NOx is formed. At the end of the inner/outer liners, the flow field becomes relatively uniform. In the 
w-velocity profiles (Figs. 13 (g) - (h)), a strong swirling flow around the fuel nozzles is observed. It is found 
that the strength of the swirling flow is slightly stronger in the main dome than in the pilot dome, which 
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Table 4. Flow area distribution!” (* Note for the cold flow calculation, both domes are not fueled). 


Aame Index Gas Mass flow rate [kg/s] 
Inflow Al Air 0.26 
Main dome Fl Fuel 0.00182* 
Pilot dome F2 Fuel 0.00182* 
Diffuser Bleed Sl Air - 0.018 
Pilot splash plate cooling A2 Air 0.0104 
Outer liner cooling 1 A3 Air 0.0053 
Outer liner cooling 2 A4 Air 0.0053 
Outer liner trim cooling A5 Air 0.0018 
Outer liner cooling 3 A6 Air 0.0024 
Main splash plate cooling A7 Air 0.0116 
Inner liner cooling 1 A8 Air 0.0096 
Inner liner cooling 2 A9 Air 0.0056 
Inner liner trim cooling Al0 Air 0.0018 
Outer liner cooling 3 All Air 0.0024 
Centerbody outer cooling | A12 = Air 0.0018 
Centerbody mid cooling Al3 Air 0.0024 
Centerbody Inner cooling | Al4 = Air 0.0024 


should be related to the fact that the airflow coming from the diffuser is aligned to the main dome (i.e., the 
pressure in the outer liner is higher). Also, the AUSM scheme predicts a higher magnitude of w-velocity 
(stronger swirling flow) than the JST scheme does. 

Figures 14 show the Mach number contour (blue: Mach = 0 and green: Mach = 0.2) at two different axial 
locations (x=0.175 [m] and 0.18 [m] ) just downstream of the pilot dome. We observe that the swirling flow 
attenuates more quickly for the JST scheme due to the numerically dissipative nature. The same conclusion 
is drawn from the Mach number contour near the main dome (not shown). 

The three-dimensional features (iso-surface of the u-velocity (green: -30 [m/s]) and v-velocity (red:-100 
[m/s], and blue: 100 [m/s]) of the predicted PVC and the dilution airflow are shown in Figs. 15 (a) JST 
scheme and (b) AUSM scheme. For the AUSM scheme, the CRZ (green) is much larger and extends upto 
the location where the dilution airflows meet. In contrast, the JST scheme predicts a much smaller CRZ, 
and the CRZ and the dilution airflows are weakly interacted. This is consistent with what we see in the 
w-velocity contour in Figs. 13 (g) - (h). 

In order to investigate the CRZ and PVC, we perform the quantitative comparison of u-velocity and 
pressure profiles in the vicinity of the main dome and pilot dome. Figs. 16 (a) and (b) show the u-velocity 
profiles in the axial direction (please see the inserted figure). We compare the location where the profile 
reaches u-velocity = -30 [m/s] and observe that the AUSM scheme (black circle) shows a much larger CRZ 
than what the JST scheme (red circle) does. More importantly, the maximum value of the reverse flow 
obtained by the AUSM scheme is -60 [m/s] (main) and -53 [m/s] (pilot), which are much stronger than the 
ones by the JST scheme (-40 [m/s] (main) and -38 [m/s] (pilot)). Under the same numerical setting (e.g., 
mesh, boundary conditions, etc.), this difference is significant. Figures 17 (a) and (b) show the gage pressure 
profiles just next to the main and pilot domes. Especially for the main dome region, the AUSM scheme 
predicts much smaller pressure, which indicates a much stronger PVC. 

Based on these observations, we can expect that the JST scheme and the AUSM scheme would predict 
flame structures in a different manner, which is our ongoing work. 
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Figure 13. Contours of gage pressure, u-velocity, v-velocity and w-velocity at the mid-plane using the JST 
scheme (left) and the AUSM scheme (right). 


V. Conclusion 


In this work, we have performed detailed validation studies of the newly implemented numerical schemes 
and applied them to the real combustor configuration, E? combustor. Compared with the available analytical 
solution & experimental data, the enhanced OpenNCC shows satisfactory performance. In the study of the E% 
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Figure 14. Contours of Mach number at x=0.175 [m] and 0.18 [m] (blue: Mach = 0 and green: Mach =0.2). 
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(a) | (b) 


Figure 15. Iso-surface of u-velocity (green: -30 [m/s]) and v-velocity (red:-100 [m/s], and blue: 100 [m/s]) 
using (a) JST scheme and (b) AUSM scheme. 
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Figure 16. U-velocity profile near the main dome (a) and the pilot dome (b) 


combustor, which is one of the most complicated CFD geometries publicly available, the preliminary results 
show that the choice of the numerical scheme might play a critical role in the accurate flow predictions, 
which could significantly influence the flame structure in a reacting flow case. From our exercise, we observe 
that the AUSM scheme is less dissipative and seems to capture the swirling flow and the recirculation buddle 
more realistically than the JST scheme. However, in terms of computational cost, the JST scheme is less 
expensive using a larger time step without numerical instability. Future work will involve the sensitivity 
analysis of combustion characteristics on the numerics and mesh quality in the E? combustor. 
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